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Abstract. We study the two-dimensional Edwards- Anderson spin- glass model using 
a parallel tempering Monte Carlo algorithm. The ground-state energy and entropy 
are calculated for different bond distributions. In particular, the entropy is obtained 
by using a thermodynamic integration technique and an appropriate reference state, 
which is determined with the method of high-temperature expansion. This strategy 
provide accurate values of this quantity for finite-size lattices. By extrapolating to the 
thermodynamic limit, the ground-state energy and entropy of the different versions of 
the spin-glass model are determined. 
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1. Introduction 

The study of disordered and frustrated systems is a current subject in statistical 
mechanics, and spin-glass models occupy a privileged place [H 121 El H]. At low 
temperatures and even in the ground state (GS), these complex systems display the 
main characteristics which dominate its physical behavior. GS quantities such as the 
domain-walls energy, for example, are frequently calculated to determine if a finite 
critical temperature exists [HI El 13 El El ITO] - 

The energy and entropy are other GS observables that give valuable information 
of these systems. The former can be calculated using different optimization techniques 
p!T| [T2]. among others, genetic algorithms [13], simulated annealing ^14j, multicanonical 
ensemble [15] and parallel tempering [161 113 IHI 119] • On the other hand, as entropy 
calculations are more difficult to carry out, sophisticated algorithms have been designed 
to determine this quantity. For example, transfer matrix [201 EH El] and ballistic- 
search [23] methods. A more popular technique to calculate entropy is known as 
the thermodynamic integration method [2H ESI E6]. This relies upon integration of 
the internal energy as function of temperature along a reversible path. Initial point 
corresponds to an arbitrary but known reference state, while the final point corresponds 
to the state for which the entropy value is required. In the practice, a problem to 
implement this technique is the necessity to calculate a suitable reference entropy at a 
not very high temperature. 

In this work we use a parallel tempering algorithm to determine both the GS energy 
and, by means of the thermodynamic integration method, the GS entropy of the two- 
dimensional (2D) Edwards- Anderson spin-glass model [27], a paradigmatic disordered 
and frustrated system. By considering different bond distributions, the more efficient 
strategies to implement the thermodynamic integration method are determined in each 
case. These strategies are based on the construction of reference states by the method 
of high-temperature expansion. We stress that, the main objective of this work, is to 
show that this implementation of the thermodynamic integration technique permits to 
obtain reliable values of the GS entropy for different disordered and frustrated models. 

The article is organized as follows. In Section El the Edwards- Anderson spin- 
glass model and the parallel tempering algorithm are presented. The thermodynamic 
integration method is described in Section El and also the construction of reference states 
is discussed. In Section IH we present our main results and the values of GS energy 
and entropy obtained by extrapolating to the thermodynamic limit from finite-lattice 
calculations. Conclusions are summarized in Section El 
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2. Model and algorithm 

We start by considering the Hamiltonian of the d- dimensional Edwards- Anderson spin- 
glass model [27] . 

H = - ^Jijaiaj, (1) 

where the sum runs over the nearest neighbors of a hypercubic lattice of linear dimension 
L and the variables (Xj = ±1 represent to N = L'^ Ising spins. The coupling constants 
or bonds, Jij, are independent random variables chosen from a given probability 
distribution. Along this work, all 2D samples (particular realizations of bond disorder) 
with a square geometry were generated with periodic-periodic boundary conditions. 

We study the Edwards- Anderson model with different bond distributions. On one 
hand, we consider the typical continuous Gaussian, 

^G(J.,) = ^exp(-4/2), (2) 
and discrete bimodal, 

MJ^,) = \ - 1) + KJr3 + 1)] , (3) 

distributions, for which the mean value is zero and the variance is one. These are the 
most popular bond distributions. To avoid confusions, we will denominate EAG and 
EAB to the versions of the Edwards- Anderson model where interactions are drawn from, 
respectively, equations ([2]) and ([3]). 

On the other hand, we also study a continuous uniform distribution with zero-mean 
value and variance one, 

^"^^^^■^"\0 if|J.,|>v^, 
and a 2p-delta distribution, 

MJ^3) = ^ E - ^b) + 6{J,, + zb)] , (5) 
i=i 

which is a generalization of ([3]). Here, b{p) is a function of p chosen so that the mean 
value and the variance of Jij to be zero and one, respectively. By direct integration it 
is easy to show that 



For p = 1 both distributions (jS]) and (jS]) are equivalent. On the other hand, 6max = 
pb^ y/3 (the maximum possible value of Jij) when p — t- oo and the equation ([5]) tend 
to the uniform distribution (jij). 

In addition, an asymmetric distribution |28] . 

PAiJ^,) = lm,-l)+5iJ^,+a)], (7) 
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where < a < 1 is a rational number, and an irrational distribution |29j . 




where c = (v^ — l)/2 ^ 0.618 is the golden ratio conjugate (or silver ratio), were also 



considered. As before, we will denominate EAU, EAP, EAA and EAI to, respectively, 
each one of these versions of the Edwards- Anderson model. 

In order to simulate these models, we use a parallel tempering algorithm [161 E] • 
It consists in making an ensemble of m replicas of the system, each of which is at 
temperature (Ti > Tk > Tm)- The basic idea of this algorithm is to simulate 
independently each replica with a standard Monte Carlo dynamics, and to swap 
periodically the configurations of two randomly chosen replicas. The purpose of this 
swap is to try to avoid that replicas at low temperatures get stuck in local minima. 
Thus the highest temperature, Ti, is set in the high-temperature phase where relaxation 
time is expected to be very short and there exists only one minimum in the free energy 
landscape. The lowest temperature, T^, is set in the low-temperature phase. 

In order to implement this algorithm, we choose equally spaced temperatures, i.e. 
Tfc— Tfc+i = (Ti — Tm) I (m— 1). Each replica is independently simulated by a single spin- 
flip dynamics where updates are attempted with a probability given by the Metropolis 
rule [30]. On the other hand, the trial exchange of two configurations Xk and Xk' 
(corresponding to the A;-th and fc'-th replicas) is attempted and accepted with probability 



where A = (/3fc/ — /Jjt) [H{Xk) — H{Xki)] and = l/T^ (without loss of generahty, we 
take the Boltzmann's constant equal to one). A unit of time or parallel tempering step 
(PTS), consists of a number ofmxN elementary spin-flip attempts followed by only one 
swap attempt. As in Ref. [17], we restrict the replica exchange to the case k' = k -\- 1. 

The parallel tempering algorithm can be used as a heuristic to obtain GS 
configurations [T8t [T9] . For this application it is not necessary to reach equilibrium, 
because only low-energy configurations are sought. Then, we have used this algorithm 
to calculate the GS energy. As in reference ^121, where only EAB and EAG models 
were studied, we have chosen m = 20 and the extreme temperatures as Ti = 1.6 and 
Tm = 0.1. In addition, the number of PTSs used here for the models with discrete 
(continuous) bond distributions, are the same ones that were used in reference [19] to 
calculate the GS energy of the EAB (EAG) model. 

On the other hand, to calculate the GS entropy we have proceeded differently. 
After an appropriate number of PTS, the parallel tempering algorithm allows to reach 
equilibrium and to calculate the mean energy at all temperatures. As we will see in 
the following section, integrating this curve one can obtain a reliable value of the GS 
entropy. 



m 




for A < 
for A > 0, 



(9) 
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3. Thermodynamic integration method 

In the following, we briefly describe the thermodynamic integration method (TIM) 
[211 [251 [26]. Given a model with a fixed number of entities (spins) N, we can write 
the basic relationship 

1 rds^ 



where U{N, T) and S{N, T) are the total internal energy and entropy, respectively. 
Integrating this equation, the GS entropy of a sample x is 

s.{N,Q) = s^{N,Tr)+ / — ^. (11) 

Ju^{N,Tn) J- 

Here, lower case letters denote quantities per spin and Tr is the temperature of a 
reference state for which entropy is known. While the integral in the second term can 
be accurately estimated by a Monte Carlo simulation, appropriate reference states are 
difficult to find. However, for an Ising spin system when Tr — > oo, a trivial reference 
state with entropy 

In 2^ 

hm s,(iV,T«) = — -=ln2, (12) 

Tn^oo iV 

is valid for any sample. In practice, good results are obtained with equation (IT^ if the 
energy is calculated for a great number of high temperatures [26] . 

Nevertheless, the TIM performance can be improved if a suitable reference entropy 
is determined at a not very high temperature. Such calculation can be made by the 
high-temperature expansion method [H [311 [32]. Let us consider the standard identity 
for Ising spin systems 

exp {f3Jij(Ti(Tj) = cosh {(iJij) [1 + djCij tanh (/? Jjj)] . (13) 

Then, the partition function for a particular sample can be write as 

= ^ JJ exp {pjijaiaj) 

= Y\ cosh {f3Jij) [1 + (TiCTj tanh {(3Jij)] 

= 2^ n ^osh (^j^^) ^ E n 11 + ^^^^ t^^h (z^^^^)] ' (14) 

where the sum run over the 2^ configurations of the system. Then, the free energy per 
spin, /, is 

-/3/ = ^[ln^xL 



ln2 + — 

N 



In j Y\ cosh {/3Ji. 



+ 



Ground-state energy and entropy of the two-dimensional Edwards-Anderson spin-glass model with different I 



+ 



N 



1^ E n [1 + ^^^^ ta^h (z^-^^^)] I ' (15) 

where [■ ■ ■]^^ represent a disorder average. Given a probability bond distribution P{Jij) 
and considering that hypercubic samples have dN bonds, the second term in the right- 
hand side of equation ( |T5l) can be write as 

/.I 

N 



In JJcosh(/3J, 

In [cosh (/3 J,,)] P{J,j)dJ, 

iij) — 

/oo 
In [cosh (iSJij)] P{Jij)dJij. 
-oo 



(16) 



Neglecting the third term in the right-hand side of equation (fTSjl . the free energy can 
be approximated by 

-/3/^ln2 + /. (17) 

From this equation, we can calculate the internal energy 

and entropy 



(19) 



Then, the integral (fT6!) should be calculated for each model to determine the entropy of 
a suitable reference state. 

For the EAB model the integral f|T6|) is 

Jb = rfln[cosh/3] , (20) 

and 

Mg^— (itanh/3 (21) 

SB ~ ln2 + dln[cosh/3] -ci/3tanh/5. (22) 

It is easy to generalize this result for the case of the discrete bond distributions such as 
the 2p-delta (E]), 



d 

Mp ~ — > ifetanh (ibl3) 
^ i=i 

d ^ d[5 ^ 
sp In 2 + - ^ In [cosh {ihp)] ^ ih tanh (i6/3) 

^ 1=1 ^ i=i 



(23) 
(24) 
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the asymmetric ([7j), 

Ma [tanh(/3) + atanh(a/3)] (25) 

SA ~ In 2 + ^ In [cosh cosh (a/3)] - ^ [tanh (/3) + a tanh (a/3)] , (26) 
and the irrational dHl), 



ui^-- [tanh (/3) + ctanh (c/3)] (27) 

si ^ In 2 + ^ In [cosh (/3) cosh (c/3)] - ^ [tanh (/3) + ctanh (c/3)] . (28) 

Notice that the two last models have the same reference state. 

For continuous bond distributions we have used the following Taylor expansion 



In [cosh (y)] = iy^- 1,^ + 1/-^ 



2/' + 



To obtain a reference state for the EAG model, let us consider the integral 

1 



FJa) 



'2ti J 



exp(-ay )y''dy 



(29) 



(30) 



It is easy to show [33j that, for integer n > 0, -fo(Q;) = l/v2a, Fi{a) = and the 
following recurrence relation holds 



'2ti J- 



exp(-ay^)y" "^dy 



dFn-2{a) 
da 



(31) 



Using the Taylor expansion (|29|) and the previous relation, we can approximate the 
integral (IT^ by the first four terms 



d 



T ^ 3 M 



and the energy and the entropy can be expressed as 

17 



mg ~ —d 



and 



sg ~ In 2 — ci 



2^ 4^ 3^ 



119 
~2A 



(3' 



(32) 



(33) 



(34) 



respectively. On the other hand, considering the Taylor expansion ( 129|) and by a direct 
integration, for the continuous uniform distribution we have that 



d 



l/32_ A/34 + A/36 _iIo8^ 

2^ 20' 35^ 280^ 



-d 



' 5' 25' 35^ 



and 



su ~ In 2 — d 



l/32_ A/34 + 3 6_ 17 8 
2' 20^ 7^ 40^ 



(35) 
(36) 

(37) 
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Previous reference states should be used carefully to calculate the GS entropy. The 
main problem arises when the samples are generated plj. Let us consider for example 
the EAB model. A canonical approach is implemented when samples are built with 
half of the bonds of each sign. We refer to these as canonical samples. On the other 
hand, if a grand canonical approach is used, grand canonical samples are generated 
in which bonds are put on each edge of the lattice with a probability given by the 
corresponding distribution (in this case ±1 bonds are chosen with equal probability). 
Figure [U^a), shows for the 2D EAB model, a comparison between the mean energy 
curves calculated with parallel tempering, for a canonical (#1) and a grand canonical 
(#2) samples (simulation parameters are given below). Here 



is the mean energy, where (■ ■ ■)t represents a thermal average. Although the lattice 
size L = 8 is small, above T ^ 4 both curves match very well with the equation f l2ip . 
Then, whatever the canonical or the grand canonical approach be used, an accurate 
GS entropy will be obtained using the reference state given by equation ( l22l) at this 
temperature. This is due to that all bonds have the same magnitude and the hyperbolic 
cosine is a even function. Therefore the integral ( IT6l) is not sensitive to fluctuations in 
the bond's values. 

A different situation arises for the systems for which only grand canonical samples 
can be generated. Figure [U^b) shows the mean energy for three distinct samples of the 
EAG model of size L = 8. Notice that the equation ( 133|) matches very well with the 
curve corresponding to sample #1, but not with those of samples #2 and #3. Then, 
the entropy only can be used as a reference state for sample #1. The problem is 
that, for this model, it is not possible to calculate a general low-temperature reference, 
because the integral ( fT6|) is very sensitive to particular realizations of disorder. Although 
these discrepancies disappear for bigger L, such samples are difficult to equilibrate at low 
temperatures. Then, to extrapolate to the thermodynamic limit, we need to calculate 
the entropy for small lattice sizes. 

A simple solution consists in to calculate a particular reference state for each 
sample. Instead of performing a disorder average, we consider the complete free-energy 
expression 




where again we have neglected the third term in the right-hand side of this equation. 
Then, the mean energy and entropy equation f|T8l) and f|T9l) . respectively, are 




(40) 



and 




(41) 
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Figure 1. The mean energy as function of T for different 2D samples of L = 8 and 
the corresponding analytical approximations (see text), (a) The EAB model and (b) 
the EAG model. 



Given a specific sample, the sums in equations (l40l) and (14T1) can be evaluated directly. 
Figure [TJ^b) shows the function (HOl) for the three samples of the 2D EAG model. From 
this comparison, it is evident that equation fj4Tl) will be a good entropy reference state 
for these samples. In fact, equations (HOj) and (HTl) can be used for all, canonical and 
grand canonical samples. 

Although both approaches can be used for discrete bond distribution, in this work 
we have only studied canonical samples of the EAB and EAA models. On the other 
hand, grand canonical samples were considered for the EAP, EAI and the models with 
continuous bond distribution, EAG and EAU. In each case, to determine an appropriate 
reference temperature Tji, first we have calculated the exact GS entropy for several 
samples up to L = 8, using a branch-and-bound algorithm [35]. Then, for these same 
samples the TIM, improved with a suitable reference as we discussed before, was used 
to calculate numerically each one of the corresponding GS entropies. Tr was chosen as 
the minimum temperature at which, for each one of the samples with 3 < L < 8, the 
numerical estimation of the GS entropy agrees, within the simulation error, with the 
exact value. 

We calculate the integral (fTTl) between this reference temperature and Tq, a very 
low temperature close to T = (the lowest temperature in the parallel tempering should 
be close to but not exactly equal to zero). Tabled] shows, for each one of the models, the 
most important parameters that we have used in the our simulations: the temperatures 
Tr and Tg, and the maximum lattice size studied, L^ax- In all cases we have used 
m = 300 replicas, 10^ PTSs (for samples with L = Lmax) and we have chosen the 
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extreme temperatures in the parallel tempering as Ti = Tr and = Tq- For lattice 
sizes L < Ljnax, an equal or smaller number of PTSs were used. 



Model 


Tr 


Tg 




EAB 


7 


0.1 


20 


BAG 


5 


0.01 


20 


EAU 


6 


0.01 


14 


EAP 


7 


0.1 


14 


EAA 


5 


0.05 


14 


EAI 


9 


0.05 


14 



Table 1. Reference temperatures used to calculate the GS entropy and maximum size 
studied for each model. 



4. Numerical results 

In this section we present the main results of our simulations. For each model, the 
disorder average for the maximum lattice size was performed over 10^ samples while, 
for smaller sizes, up to 10^ samples were necessary to obtain accurate values of the GS 
energy and entropy. 

To extrapolate, we have fitted our data of the GS energy to this scaling function 

fu = Uoo + guL-'^^ (42) 

where u^o is the GS energy value in the thermodynamic limit and g^ and du are two 
parameters. A similar scaling function 

fs = s^ + gsL-''% (43) 

with parameters gs and (i^, was used to estimate the thermodynamic limit of the GS 
entropy, Soo- For a given fit, the corresponding goodness-of-fit parameter Q is calculated 
[36] . A value Q > 0.1 is considered as indication of good quality of the fit. As is usually 
the case, scaling functions do not include all possible finite size corrections, and therefore 
better fits are obtained when data for very small sizes are left out. On the other hand, 
leaving out too many points can result in large error bars for the best fit parameters. 
Then, the results presented here were obtained by fitting the data over the largest range 
that gives a goodness-of-fit of Q > 0.1. 

4.1. EAB and EAG models 

Figure [2] (a) shows the GS energy for the EAB and EAG models and different lattice 
sizes. The results of fitting these curves with the scaling function fH2|) are presented 
in table El Although (for these models) the maximum lattice size simulated are not 
very large, the values of Uoo agree very well with the most accurate results reported in 
previous works: u^o = -1.40193(2) [37|, Moo = -1.40197(2) [38] and u^o = -1.4009(3) 
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Figure 2. The GS (a) energy and (b) entropy versus L, for the EAB and EAG models. 
The dotted lines are the best fits obtained. The insets show, for both models, the GS 
energy and entropy as function of N^^. 



[H] for the EAB model; = -1.31479(2) [38] and = -1.3149(5) [19\ for the 
EAG modeL Only the exponent du = 1.95(11) for EAB model is a little different from 
du = 2.13(4), the value that we have obtained previously [19j. The reason for this 
discrepancy is that the fit of the GS energy with equation fH2|) . is very sensitive to the 
range of L used (in reference |I9] was used a range of 5 — 30). 

The GS entropy for these same models is shown in figure [2] (b), while in table [3] are 
presented the results of fitting these curves with the scaling function ( H3l) . The value 
of Soo = 0.0714(9) for the EAB model is close to the results reported in the literature: 
= 0.075(2) [39], Soo = 0.0709(6) [lO] and s^ = 0.078(5) [23]. On the other hand, 
because the EAG model has a non-degenerated fundamental level (that has only two 
configurations related by a global spin-flip), the thermodynamic limit of the GS entropy 
is expected to be zero. However, we obtain Soo = 0.0003(1), a value very close to, 
but not zero. This is due to the existence of systematic errors in the implementation 
of the thermodynamic integration method [26]. Among others, one source of error is 
the (bad) assumption of that the temperature Tq is sufficiently low, so as not to affect 
the calculation of the integral (ITT]) . To improve our result, we should equilibrate each 
sample up to a lower temperature. As this is very hard to do, we assume that the 
entropy values are affected by a systematic error of order 10^. This shows the accuracy 
of our implementation of the thermodynamic integration method, and justifies why a 
system with a trivial GS as the EAG model was studied. Finally, in table |3] we can see 
that the exponent dg for this model is very close to 2. This is correct because the exact 
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Figure 3. The GS (a) energy and (b) entropy versus L, for the EAU model and the 
EAP model with different values of p. The dotted lines are the best fits obtained. The 
insets show the GS energy Uoo and entropy Soo as function oi p~^. 



GS entropy is simply 

s = 4 In 2. (44) 

The inset in figure [2] (b) shows that this equation and the GS entropy for the EAG 
model agree very well for all sizes. 

4.2. EAU and EAP models 

Figured (a) shows the GS energy for the EAU and EAP models for different lattice sizes 
and p values. The inset shows the dependence with the parameter p of the GS energy 
Uoo, which is obtained by fitting the data with the scaling function fH2]) . We can see 
that for p > 2, the curve monotonically tends to the corresponding value for the EAU 
model (see table [2]). The same is observed in figure |3](b) for the GS entropy. 

By fitting the entropy curve of the EAU model, we obtain a negative value 
Soo = —0.0003(1). As this model has a continuous bond distribution (jl]), it should 
have a non-degenerated fundamental level with entropy zero. Again, the problem is 
the systematic error affecting the results produced by the thermodynamic integration 
method. Nevertheless, figure [3] (b) shows that the GS entropy calculated with this 
technique for each lattice size follows the exact curve (jH]). 
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4.3. EAA model 

In the EAA model, bonds are chosen randomly with equal probability between two 
values: +1 (ferromagnetic bond) and —a (antiferromagnetic bond). We could have 
defined a similar distribution by exchanging the signs of these bonds, 

PA*{-hi) = \ W^, - a) + (5(Ji,- + 1)] . (45) 

We will call this the EAA* model. Figures H] (a) and (b) show, respectively, the GS 
energy and entropy for the EAA and EAA* models for two different values of parameter 
a. The most outstanding behavior is that both quantities seem to oscillate for the EAA* 
model. In particular, the GS energy and entropy for the two models, do not match for 
odd values of the lattice size. This discrepancy increases with decreasing a, but tends 
to disappear with increasing L. 

The reason of this behavior is in the geometry of the lattice, which is square, and 
the periodic-periodic boundary conditions. For small values of a, antiferromagnetic 
bonds in the EAA model are weak and the GS is formed by ferromagnetic islands. This 
structure is not affected considerably by the L value: ferromagnetic clusters can fill the 
lattice with an even or odd L. However, for the EAA* model the GS is antiferromagnetic 
and a different situation arises. Now, due that nearest-neighbors spins point in opposite 
directions, the GS configurations are characterized by the homogeneous magnetization 
of both sublattices but with different orientations. The GS energy and entropy are 
affected by the periodic-periodic boundary conditions, because to accommodate a large 
(percolating) antiferromagnetic island in a lattice with odd L, it is necessary to create 
an energy wall of the same length order. This does not happen for lattices with even L, 
for which the antiferromagnetic structure can fill the system. As the defect energy (or 
entropy) depends approximately on L and the GS energy (or entropy) depends on L^, 
the discrepancies disappear with increasing size and, in the thermodynamic limit, both 
models will have the same u^o (soo)- 

In order to carry out the fitting, we have only used the EAA data. Insets in figures 
m (a) and (b) show, respectively, the dependence with parameter a of the GS energy 
Moo and entropy s^. For a = 1, the EAA and EAB models are equivalent and therefore 
have the same values of energy and entropy. When a = half of the bonds are zero and 
the other half —1, and then the GS energy is lima^o^oo = —1- The GS entropy is also 
easy to calculate: since the probability that one spin has four bonds of zero strength is 
(1/2)^ = 1/16, in the thermodynamic limit the number of free spins (those spins whose 
fiipping does not change the energy of the state) will be n = A^/16. Then 

lim = 4 In (2") = ~ 0-0433. (46) 

a-S-O N 16 

The insets in figures H] (a) and (b) show that the GS energy changes smoothly 
between these limits but the GS entropy does not. The behavior of Soo can be explained 
considering that the number or free spins depends strongly on the parameter a. For 
a = 1 all bonds have the same magnitude (i.e. ±1) and spins with two frustrated and 
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Figure 4. The GS (a) energy and (b) entropy versus L, for the EAA and EAA* 
models, for two values of a as indicated. The dotted lines are the best fits obtained. 
The insets show the GS energy Uoo and entropy Soo as function of the parameter a. 



two satisfied bonds will be free. If the parameter diminishes a little, a < 1, n changes a 
lot: in order to have a free spin now it is necessary that, either four bonds converging 
onto it have the same magnitude (with two of them being frustrated and two satisfied), 
or that two of them have magnitude 1 and the other two a (with two bonds of magnitude 
1 and a frustrated and the other two bonds satisfied). Because the probability of this 
happening is smaller than before, the GS entropy falls abruptly [see inset in figure |l](b)]. 
A new possibility arises when a = 1/3: in a GS configuration, a spin with three satisfied 
(frustrated) bonds of magnitude a and one frustrated (satisfied) bond of magnitude 1, 
will be free. Then, the GS entropy should increase a little at a = 1/3. 

Another characteristic of the EAA model, is how the energy gap AHq between the 
GS and the lowest excitation state depends on the parameter a. Supposing that these 
excitations are due to single-spin flips only, at a = 1 the gap is AHq = 4 and this 
corresponds to flipping a spin with three bonds satisfied and one frustrated. Also, for 
1/3 < a < 1 is AHq = 2(1 — a) [this gap correspond, for example, to spins with three 
bonds of magnitude a, with two of them being frustrated and one satisfied, and the 
remaining bond of magnitude 1 satisfied], for l/5<a<l/3is AHq = 2(1 — 3a) [spins 
with three frustrated bonds of magnitude a and one satisfied bond of magnitude 1] and 
for a < 1/5 is AHq = 4a [spins with four bonds of magnitude a, with three of them 
being satisfied and one frustrated]. Interestingly, if we consider excitations in which are 
involved many spins, for a < 1 there are a few samples with a smaller gap. For example, 
at a = 0.1 the gap should be AHq = 0.4, but we have found samples of size L = 10 with 
AHq = 0.2: excitations are due to droplets whose walls are formed by a net number of 11 
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Figure 5. The GS (a) energy and (b) entropy versus L for the EAI model. The dotted 
lines are the best fits obtained. The inset in (a) shows the histogram of the energy gap 
obtained for 10^ samples oi L = 10. The inset in (b) shows a comparison between the 
GS entropy s^o of the EAI, EAB and EAG models, as function of A^~-^. 



satisfied bonds of magnitude a and one frustrated bond of magnitude 1. Nevertheless, 
as the parameter a is a rational number, the energy levels of this model form a discrete 
spectrum. 

4.4. EAI model 

Figures O (a) and (b) show the GS energy and entropy for the EAI model and different 
values of L. The results of fitting these curves with the scaling functions P2|) and (143!) 
are presented, respectively, in tables [2] and El 

Due that the parameter c in the bond distribution (jH]) is an irrational number, it 
is expected that the energy levels form a dense spectrum [41j. We have found evidence 
that indicates that this is correct. Inset in figure [5] (a) shows a histogram of the energy 
gap for the EAI model obtained for 10^ samples of L = 10. Although most of the 
samples have a AHq ^ 0.76393 (corresponding to excitations of a single-spin flip), a 
considerable number of them have a smaller energy gap (corresponding to big droplets). 

Finally, we can see that contrary to the EAG and EAU models, which have also a 
dense spectrum of energy levels, here the GS entropy is not zero. Inset in figure [5] (b) 
shows a comparison between the curves of Soo versus A^~^, for the EAI, EAB and EAG 
models. It is evident that the GS of the EAI model is degenerated. An exponential 
number of GS configurations exists because the bond distribution Pi is discrete. As 
the magnitude of bonds are similar to those of the Pa distribution with a = 0.618, the 
number of free spins and the GS entropy for both models should be similar. We obtain 
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Model 




Qu 


du 


Range of L 


EAB 


-1.4024(10) 


1.2(2) 


1.95(11) 


5-20 


EAG 


-1.3136(13) 


1.5(4) 


2.39(22) 


4-20 


EAU 


-1.3763(8) 


1.6(1) 


2.32(6) 


3-20 


EAI 


-1.1713(6) 


1.6(2) 


2.78(7) 


3-20 



Table 2. Best fit parameters for the scaling function p2)) and tire range of L used. 



Model 




9s 


ds 


Range of L 


EAB 


0.0714(9) 


0.71(7) 


1.56(6) 


5-20 


EAG 


0.0003(1) 


0.694(7) 


2.00(1) 


3-20 


EAU 


-0.0003(1) 


0.70(8) 


2.01(1) 


3-14 


EAI 


0.0209(5) 


1.01(4) 


2.12(3) 


3-14 



Table 3. Best fit parameters for the scaling function (j43|) and the range of L used. 

Soo = 0.0209(5) for the EAI model, which are very close to s^o ~ 0.017 for the EAA 
model (with a = 0.618). 

5. Conclusions 

In this work we have used a parallel tempering algorithm to determine both the GS 
energy and, by means of the thermodynamic integration method, the GS entropy of the 
2D Edwards-Anderson spin-glass model with different bond distributions. 

To implement the thermodynamic integration technique, we have built reference 
states by the method of high-temperature expansion. Although different strategies can 
be used for canonical and grand canonical samples a simple solution, consisting in to 
calculate a particular reference state with equations (HO!) and (14T!) for each sample, 
works in all the cases. By using this method, we have been able to calculate accurate 
values of the GS entropy. This allowed us to make an study of six versions of the 2D 
Edwards-Anderson spin-glass model, which have different GS properties. 
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